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Because of the nonlinearity of the Einstein equations, the cosmological fluctuations which are 
generated during inflation on a wide range of wavelengths do not evolve independently. In particular, 
to second order in perturbation theory, the first order fluctuations back-react both on the background 
geometry and on the perturbations themselves. In this paper, the gravitational back-reaction of long 
wavelength (super-Hubble) scalar metric fluctuations on the perturbations themselves is investigated 
for a large class of inflationary models. Specifically, the equations describing the evolution of long 
wavelength cosmological metric and matter perturbations in an inflationary universe are solved to 
second order in both the amplitude of the perturbations and in the slow roll expansion parameter. 
Assuming that the linear fluctuations have random phases, we show that the fractional correction 
to the power spectrum due to the leading infrared back-reaction terms does not change the shape of 
the spectrum. The amplitude of the effect is suppressed by the product of the inflationary slow-roll 
parameter and the amplitude of the linear power spectrum. The non-gaussianity of the spectrum 
' induced by back-reaction is commented upon. 
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I. INTRODUCTION 



The study of cosmological fluctuations is one of the cornerstones of modern cosmology. In order for a cosmological 
model to be considered successful, it must be able to reproduce, among other things, the power spectrum of the 
perturbations. These perturbations leave their mark as anisotropics in the cosmic microwave background (CMB) and 
ly-j ■ go on to act as seeds for structure formation. The theory of cosmological perturbations establishes the bridge between 
observations (namely observations of fluctuations in the CMB and in the distribution of structure in the universe) and 
the physics of the very early universe which is responsible for providing the generation mechanism for the fluctuations. 

Arguably, the most successful cosmological paradigm for the early universe is the inflationary scenario . Inflation- 
ary cosmology provides a mechanism for generating the primordial fluctuations and predicts a nearly scale-invariant 
Q H l spectrum of Gaussian, adiabatic, scalar fluctuations Q (see also These predictions have been verified ex- 

perimentally to high accuracy in the last decade (see, e.g., the most recent WMAP CMB anisotropy maps |||). Since, 
in addition, inflationary cosmology solves numerous problems of Standard Big Bang cosmology [lj (for a review of 
these problems along with a discussion of some conceptual problems inherent to inflation see 0), inflation is widely 
C$ . considered as a key ingredient of the theory of the early universe. 

At the present time, however, inflationary cosmology does not quite have the status of a theory. It is best thought 
of as a successful scenario that resolves many of the problems that plague Big Bang cosmology. There are a large 
number of different models that result in accelerated (i.e. inflationary) expansion. However, many models share a 
common feature in that they involve a scalar field (the "inflaton" ) which undergoes a period in which it rolls slowly 
down it's potential - leading to what is known as "slow- roll inflation". Our analysis will be in the context of a general 
slow roll inflation model. 

The theory of cosmological perturbations (see e.g. for a comprehensive review, and Q for a recent abbreviated 
overview), a formalism crucial to the understanding of CMB anisotropies, is usually studied within the framework 
of linearized gravity. One writes down an ansatz for the form of the perturbed metric about a homogeneous and 
isotropic background space-time, linearizes the Einstein equations in the amplitude of the perturbations, and solves 
the resulting equations. In this scheme, all Fourier modes of the fluctuations evolve independently, as must be the 
case in a linear approximation. The Einstein equations which govern the evolution of space-time and matter are, 
however, nonlinear. Thus, retaining terms quadratic and higher in the perturbation amplitude leads to interactions 
between different perturbation modes. These interactions determine the "gravitational back-reaction" , the difference 
between the full evolution of the space-time and what would be obtained in linear theory, and will lead to potentially 
important modifications of the results obtained at linear order. In particular, they may effect the key qualitative 
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predictions of inflation, namely the scale invariance of the spectrum and its Gaussianity. 

The period of inflation in most scalar field-driven inflationary models is very long (measured in units of the Hubble 
time during inflation) . Thus, the red-shifting of scales leads to the population of a large phase space of long wavelength 
modes (modes with a wavelength larger than the Hubble radius). The back-reaction of such long wavelength modes on 
the background space-time was first studied in [1(1 1 1 11] (making use of the concept of an effective energy-momentum 
tensor of fluctuations first used in studies of short wavelength gravitational waves 12]). It was found that this effective 
energy-momentum tensor acts like a negative cosmological constant with a magnitude which increases in time as 
the phase space of long wavelength modes grows. A physical explanation of this effect in the quasi-homogeneous 
approximation to the evolution equations was recently provided in The effect can become non-perturbatively 

large if the period of inflation is sufficiently long and leads to a change in the Hubble expansion rate. This change 
is physically measurable in models with at least two matter fields |14|. In models with only one matter field, however, 
the leading infrared back-reaction terms are not physically measurable by a local observer 15. 16, 17| (sec also 18]. 
and see 191 for a review of previous work on gravitational back-reaction in inflationary cosmology [40|). 

Given that the back-reaction of first order cosmological perturbations on the background cosmology can become 
large, it is important to determine whether this gravitational back-reaction can also lead to large effects on the 
fluctuations themselves, potentially also changing their key characteristics like almost scale- invariance and Gaussianity. 

It is the intent of this article to solve the perturbed Einstein equations to quadratic order, and determine the 
modifications to the results of the linear analysis, focusing on the effects on long- wavelength fluctuations. There has 
been a significant body of previous work devoted to second order cosmological perturbations, see e.g. 
EM I24I I23 and papers quoted therein. The effect of long-wavelength modes has also been studied in the "separate 
universe" approach [l8L l26l 12^ , in which the effect of long- wavelength perturbations in encoded as a change in the 
background geometry. Our approach is similar. In particular, we will neglect spatial gradients in the equations of 
motion, and thus are focusing on the leading infrared contributions to back-reaction. What differentiates our analysis 
from previous work is the emphasis on the fact that, in an inflationary universe, modes continuously move into the 
infrared sector (wavelengths greater than the Hubble radius), and that thus the infrared phase space grows. This 
leads to the concern that back-reaction effects grow without limits, a concern which is the main motivation for our 
work. 

The main result of our study, however, is that the leading infrared contributions of back-reaction to the power 
spectrum of cosmological fluctuations is very small. Assuming that the linear fluctuations have random phases, as 
they do in the simplest inflationary models, the relative contribution of our back-reaction terms to the power spectrum 
is suppressed by the product of a dimensionless inflationary slow-roll parameter and the amplitude of the linear density 
fluctuations. 

We free ourselves of constraints imposed by any specific inflationary model by abstaining from picking a specific 
form of the inflaton potential. Rather, we assume that the inflationary slow-roll conditions are satisfied and we 
employ a very general form of the potential, one that allows us to easily interpret our results within the context of 
an explicit slow-roll realization. Thus, with regards to cosmological models, our approach is extremely general, while 
still retaining its ability to be specific. 

Section 2 of this paper is concerned with the general setup of the calculation. Here, the background equations are 
solved to second order in the slow-roll parameter (e), neglecting spatial gradients, and the corrections to both the 
background scalar and the Hubble constant are determined, again, to second order in e. Section 3 begins with a 
brief review of classical, relativistic, cosmological perturbation theory. The linear perturbation equations are solved 
using the background obtained in the previous section. The non-linear equations make their appearance in Section 4, 
where the next-to-leading order terms (those due to back-reaction) are solved for. Section 5 interprets these results as 
modifications to the linear terms, and their effect on the spectrum of perturbations is determined. Finally, the effects 
of back-reaction on the Gaussianity of the perturbations is investigated in Section 6. 

II. BACKGROUND 

Before solving for the perturbations, we must determine the appropriate background for our model. We confine our 
attention to a model consisting of pure gravity (with vanishing cosmological constant) and a single, homogeneous, 
scalar field 4>, which we presumes satisfies the slow-roll conditions. We write 

<t>{x,t) = o + e/i(t) + e 2 / 2 W, (1) 

where (f>o is a constant and e is a dimensionless slow-roll parameter whose value depends on the specific model in 
question. In consequence of our background space-time will be approximately de Sitter space, 



(t) = e ff(t)t 



(2) 
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Here, the Hubble rate H(t) is slowly time-dependent and can be expanded as 

H(t) = H + eh 1 (t) + e 2 h 2 (t), (3) 

where H a constant. 

The dynamics of the scalar field is determined by its potential V((f>). In order not to limit our results to any specific 
model, but, rather, to show them to be a generic feature of slow-roll inflation, we use the following expansion (related 
to a power series expansion of the potential about the field value </>o): 

V(4>) = A^o + eA^(t) + e 2 Acj ) 2 (t) = + + V™ , (4) 
where superscripts indicate the order in e and 

= V(ct> ) - <f>oV' (M + \<&V" (0 O ) , (5) 
eX = V' (fa) ~ <f>oV" (cj> ) (6) 
e 2 A = \v'\M, (7) 

the primes denoting a derivative with respect to the field. Thus, we sec that /i,A, and A are dimensionful parameters 
depending entirely on the form of the potential in the neighborhood of <fio- In essence, this represents the series 
expansion (in e) of any potential that can be used to generate slow-roll inflation. 
In light of this, we can use the Klein-Gordon 

+ 3-0+^=0, (8) 
a o<p 

and the Friedmann equations 

3d 2 + 47ra 2 2 + 8ira 2 V = 0, (9) 

• 2 

% + 2- + 47r^ 2 - 8nV = (10) 

(solved order by order in e) to solve for the background, accurate to second order in the slow roll parameter. 
We find that 

H 2 = 8n^ , (11) 



and 



A 1 e -3H t 

hit) = ^r(^r-*-^-)' ( 12 ) 



= T^T7^(- 2 + ™ + ^ 3Hot (3H Q t + 2)) , (13) 



h 2 (t)= - 12967rM 2 g 4 t ( 288g o 2 A 2 ^V + 162g 5 A 2 t^ - 243Aff 7 t (15) 
+16A 2 7rVe- 6Hot - 128A 2 tt Ve^ 3ff °* - 288 A 2 V 2 /i 2 H t + 112A 2 ttV) . 



Note that the corrections to the scalar field have the property that they and their first temporal derivatives vanish 
at t = 0, which coincides with the onset of slow-rolling. Despite the overall factor of ir 1 preceding h 2 {t) 1 this term is 
not singular at the origin, as can be seen by expanding in a power series about t = 0. 
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III. LINEARIZED THEORY 

The theory of cosmological perturbations is usually restricted to an analysis of the linearized (in the amplitude of 
the perturbations) Einstein equations. In this section, we review the salient results (see e.g. 0, for comprehensive 
reviews) . 

The most general perturbed line element can be written as 

ds 2 = a 2 (r))[{l + 2<J>)dr/ 2 - 2{B,, + S t )dx l dr) - [(1 - 2^) 7y - + 2E- ij + (F i;j + F j;i ) + h ,, d.r <l.r r (16) 

where ( &, V E', B, E represent scalar, Si,Fi vector, and hij tensor metric perturbations, respectively. These are distin- 
guished by their transformation properties under three dimensional rotations. For inflationary cosmology and in linear 
perturbation theory, we can discard vector modes since they do not grow, and tensor modes because they grow at a 
slower rate than scalar metric fluctuations. Hence, we focus on the scalars. It is possible to choose a gauge in which 
E = B = 0, the so-called 'longitudinal' or 'conformal Newtonian' gauge. This gauge is convenient for calculational 
purposes. For matter without anisotropic stresses to linear order in the matter field perturbations, it follows from 
the off-diagonal space-space Einstein equations that $ = 'J. Hence, for such matter the metric in longitudinal gauge 
takes the form 

ds 2 = (l + 2n<i>(x,t))dt 2 -a 2 (t)(l-2n<I>(x,t))(dx 2 +dy 2 + dz 2 ). (17) 

In the above, n is a dimensionless parameter which indicates the order of the term in gravitational perturbation 
theory. 

At second order, the perturbed metric is in general much more complicated. In particular, there is mixing between 
scalar, vector and tensor modes, and an anisotropic stress in generated, leading to $ ^ However, it can be seen 
explicitly that all of the complicating terms contain spatial gradients and can hence be neglected in our study of 
the back-reaction effects of long wavelength modes. Thus, if we work to leading order in the infrared terms, we can 
neglect vector, tensor modes and anisotropic stress. Thus, we can apply longitudinal gauge also at second order. We 
will use this gauge in the following. 

To incorporate effects due to slow-rolling, we expand $ in powers of e to obtain 

$i(x,t) = *i(f,i) + e ai (x,t) + e 2 [3 1 {x,t). (18) 

The subscript i denotes the fact that the effects are linear in n. 

In order to source our first order metric perturbations, matter perturbations 5\{x, t) must be present, i.e. 

4>{x,t) — ► (f>(t) + KSx(S,t) , (19) 

where <fi(t) is the background matter solution given by 

At the linear level, perturbations decouple in Fourier space. It is thus convenient to track the evolution of each 
Fourier mode individually. At higher orders in perturbation theory, there will be mixing between the modes. Typically, 
cosmological perturbations are classified into two distinct sectors: sub-Hubble (UV) and super-Hubble (IR) modes. 
During inflation, the phase space associated with super-Hubble modes grows exponentially, while that of the sub- 
Hubble modes remains constant. In addition, the modes which are important from the point of view of structure 
formation and CMB anisotropics are super-Hubble during the last 50 e-foldings of inflation. It is for these reasons 
that we focus our attention on these modes and ignore their UV counterparts. This choice justifies our dropping of 
spatial derivatives in the equations of motion. 

Using the standard result for the energy-momentum tensor of a scalar field 

T„„ - - g^A^p 1 - V (</>)) , (20) 

expanding the Einstein equations in a power series in k, and truncating after first order leads to the equations of 
motion for scalar cosmological perturbations, which read (see e.g. Q) 

• 2 

(ii): +8tt^'i -2$i -8-$i -4$i^r - 16tt$i((/.) 2 (21) 

a a a* 

-8ttU (1) (0) + 16ttU (0) (</>)$! = 0, 
(00): 6$i- +87nMi+87rU (1) (</>) + 167rU (0) = 0, (22) 
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where the overdot denotes a derivative with respect to cosmic time, and the {it) and (00) indicate the tensor indices 
of the Einstein equations in question. 
These equations can be solved to yield 

*i(#,t) = 0, (23) 
ai (x,t) = J d 3 fc/(fc)e^e"W 1/2 (24) 

ih@,t) = _^ 2A ^o-a 2 +e _ got) r d ^ fi zyk.3 eia%vl/ 2 (25) 

A^ J 

Sl (S,t) = -2^tl J d 3 kf(k)e^ s e^V^ 2 , (26) 

where ipi a constant representing the amplitude of the spectrum of linear fluctuations, f(k) is the function describing 
the shape of the spectrum, are the phases (assumed later on to be random), and V is the cutoff volume used in the 
definition of the Fourier transform. Note that k denotes the co-moving momentum. We have introduced the cutoff 
volume in the definition of the Fourier transform in order that the dimension of the Fourier transform <£>(fc) of $(x) 
is fc 3 / 2 , which in turn ensures that there is no volume arising in the relation between the power spectrum P$(fc) and 
$(fc) (see Section 5). 

Let us comment briefly on the interpretation of this solution. Note first that the fact that vanishes is a consistency 
check on our analysis. In the pure de Sitter limit (no rolling of the scalar field), there are no scalar metric fluctuations 
to first order in perturbation theory (in k). Since the equation of motion is second order, there are two fundamental 
solutions for each mode. In the long wavelength (super-Hubble) limit, the dominant solution is constant in time, and 
the sub-dominant solution is a decaying mode. We see that a\ is a the constant perturbation sourced by S\, while (3\ 
is generated by the rolling of the scalar field and is a combination of the constant and decaying mode. 



IV. BACK-REACTION 



To second order in gravitational perturbation theory (expansion in k), there are interactions between the Fourier 
modes of the fluctuation variables caused by the non-linearities of the Einstein equations. In order to determine the 
effects of non-linearities, we make the substitutions for the metric perturbation variable 

$(2,t) — >$i(£,t) + /e$ 2 (£,i), (27) 

where $2 can be in turn expanded in terms of the slow-roll parameter e as 

$2 0?, t) = # 2 (a , t) + ea 2 (x, t) + e 2 (3 2 (x, t) , (28) 

and for the matter field 



<p(x, t) — > <p(x, t) + k 2 S 2 (x, t) . (29) 

Note that $2 and S 2 represent the effects of gravitational back-reaction. 

In the following we will neglect spatial gradient terms in the equations of motion, since we are interested in the 
infrared modes and in the coupling between different infrared modes j^jj. To order k 2 in the perturbative expansion, 
the Einstein equations then become 

• 2 

(n): 2$2 + 8$2- + 8$2(- + 27r0 2 + ^ T -27rV r(2) (< ? ! ) ))-87r^2 (30) 
a a la 

. 2 

= 4$i($i +2$i^r +6$i- +4$i- +8tt$i(Mi + ^V (1) (0)$i, 
a a a 

(00) : 6$ 2 - + 871-^2 - 167ry (0) ((/>)$2 (31) 
a 

= 3$i 2 - 4tt5i 2 - 12$i$i- - 167i-y (1) ((/))$i - 8nV {2) U). 

a 
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These equations can be solved to yield 

Mx,t) = 0, (32) 

a 2 (x,t) = 0, (33) 

A(f,t) = ^ (A ^ a " Aa) (-l + e- g oV(a), (34) 

5 2 (x,t) = 0, (35) 

with 

g (x) = ( d 3 kf(ky ts e la *V 1/2 . (36) 



Let us briefly comment on the physical interpretation of these results. First, the vanishing of ^ 2 is a consistency 
check since there are no scalar metric fluctuations in pure de Sitter space. Since the linear fluctuations are first 
order in e, they will only contribute to the second order perturbations to quadratic order, and hence the vanishing 
of a 2 is another consistency check on the algebra. From the expression for $2 it is manifest that the second order 
perturbations are generated by the linear inhomogeneities $1 and Si at quadratic order. The vanishing of d 2 is an 
interesting and unexpected result. It says that, to this order, there are no back-reaction effects on the evolution of 
the background scalar field. 

V. EFFECTS OF BACK-REACTION ON THE POWER SPECTRUM 

Having now determined the form of the back-reaction terms, it is important to estimate their amplitude. From 
observations of CMB anisotropies, we know that the linear perturbations are of order Kipie ~ 10~ 5 . The back- 
reaction terms should be expected to be of order (tnpie) . However, the second order correction to a fixed Fourier 
mode receives contributions from all linear Fourier modes to this order. Hence, one could expect the back-reaction 
effect to be amplified by a phase space factor which measures the phase space of Fourier modes which contribute. 
Since in inflationary cosmology, the phase space of infrared modes is growing, and the linear fluctuations do not 
decrease in amplitude on scales larger than the Hubble radius, the effects of back-reaction could be expected to grow 
in time and become non-perturbatively large. In this section we show that, provided that the linear fluctuations have 
random phases, the leading infrared quadratic back-reaction effects of linear fluctuations on the power spectrum of $ 
do not show any large phase-space enhancement. 

The total power spectrum Vtotai(k) of $, including the leading infrared terms of second order in k, can be written 

as 

Ttotadk) = fc 3 |$fe| 2 = V 1 (k)+P 2 (k), (37) 

= ^(fc) + $ 2 (fc)| 2 fc 3 , 

where $i(fc) are the Fourier transforms (using the definition of Fourier transform including the cutoff volume as in 
Equation 124(1 ') of &i(x). Making use of the results of Sections 3 and 4 we have 

ix(fc) = c^i/(A;)e fa fi[l-e£(t)] 

$ 2 (fc) = e 2 ^lE{t)h{k) , (38) 
where h(k) is the Fourier transform of g 2 {x) and where we have introduced the symbol E(t) for the function 

E(t) ee 27 W° 2 - A 2 (-1 + e - Hot) m 

The leading back-reaction correction to the power spectrum, denoted by the function VBR(k), comes from the cross 
term in l(38|l and is thus linear in <& 2 (k): 

V BR {k) = 2 K |$ 1 (fc)$ 2 (fc)|fc 3 . (40) 

Thus, the fractional correction to the power spectrum due to the leading back-reaction contributions is 

V BR {k) _ e 2 tfE(t)h(k) 
= Zk- 



Vi(k) e^i/(fc) + 0(e e ) 
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from which it follows that, modulo the ratio of h(k) over f(k), the back-reaction terms are suppressed, as expected, 
by Keipx- The ratio of h(k) over f(k) is the possible large phase space enhancement factor. 

Before continuing, we specify the linear power spectrum. We choose a normalization wavenumber k n and choose 
to be the amplitude of the power spectrum at k = k n . The function f(k) describes the spectral shape. We choose 
a power law with a tilt £ away from scale-invariance, i.e. we write 

/ k \ ~ 3/2 ~ C 

m ={lt) K ■ (42) 

It can easily be checked that V\{k n ) — 

We now evaluate the magnitude of h(k), assuming that the phases are random: 

h(k) = _L_y-i/2y d 3 xg 2 (x)e- iSs 

(PkJik^fik - Z 1 )e Ka *i +a *-*i V 1/2 . (43) 

Given that we are considering the effects of long- wavelength fluctuations, we must restrict the above integral over fci 
to run only over super-Hubble modes, i.e. 

\ki\ < H . (44) 

To estimate the magnitude of h(k), we insert the spectrum (|42|l into (|43|l . If we consider the effects of back-reaction 
on modes k which are sub-Hubble now, we can apply the approximation 

in which case the integral simplifies. 

Assuming constant phases for the moment, the integral (|43[) can be easily estimated 

h(k) ~ k^/^k^H^V 1 / 2 . (46) 

Note in particular from l|4tj|) that the k-dependence of VBB.(k) is the same as that of the linear power spectrum. The 
leading effect of back-reaction thus does not change the power index of the spectrum. However, for wavelengths close 
to the Hubble radius, the approximation l|45(l is no longer good, and the correction terms will yield changes to the 
index of the power spectrum. The second fact to notice about the result (|46(l is the cutoff volume divergence. This 
stems from the fact that as V increases, more and more infrared modes are contributing to the back-reaction. For 
constant phases, the effect is additive. The volume divergence thus represents the phase space enhancement which is 
the focus of this investigation. 

Let us now consider the more realistic situation - realized in typical inflationary models - in which the phases 
are random. A simple way to estimate the effects of the random phases in (|43|) is to add up the amplitudes of the 
back-reaction contributions of all infrared modes k\ as a random walk. This means dividing the amplitude obtained 
previously by N{V) X ^ 2 , where N(V) is the number of modes. Since for a finite volume V the wavenumbers are 
quantized in units of Afc ~ V^ 1 ^ 3 , the number N(V) becomes 

N(V) ~ (A) 3 (47) 

in which case the result (|43|l becomes 

h(k) ~ k-^^k^H^ . (48) 



Inserting this into l|41|). we obtain our final result 



VBR{k) 2^/^). (49) 



Vi(k) Y \H, 

The main conclusion we draw from (|49|l is that there is no phase space enhancement of the back-reaction of long 
wavelength modes on the spectrum of cosmological perturbations, in contrast to the positive enhancement found for 
the back-reaction on the background metric. Given the absence of such a phase space enhancement, we find - as 
expected - that the back-reaction terms in the power spectrum are suppressed compared to the terms coming from 
the linear perturbations by neipi. Thus, they are completely negligible in the case of a COBE-normalized spectrum 
of almost scale-invariant linear fluctuations. In addition, we find that the leading back-reaction terms do not change 
the spectral index. 
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VI. NON-GAUSSIANITY OF THE SPECTRUM DUE TO BACK-REACTION 



Having established that the back-reaction of infrared modes cannot substantially modify the amplitude and spectral 
tilt of the power spectrum of perturbations, we make some comments regarding the effects of back-reaction on the 
Gaussianity of the spectrum. 

The inclusion of higher-order terms implies correlations between different modes, thus breaking strict Gaussianity. 
However, the question remains: how badly broken is it? To estimate this, we turn our attention to the bispectrum 
(three-point function) . 

In the case of purely Gaussian distribution, all odd moments are identically zero. Therefore, the non-vanishing of 
the bispectrum indicates that the distribution cannot be Gaussian. 

We take it for granted that the bispectrum does not vanish (for examples of the three-point function see [2fil l2Sj . 
In the context of higher order perturbation theory (see e.g. |21|), however, its amplitude is exceedingly small. We 
estimate it to be no larger than of order K 4 e, thus making it quite unlikely to be detected experimentally. Thus, we 
conclude that, although back-reaction modifies the distribution, Gaussianity remains an excellent approximation. 



VII. CONCLUSIONS 



In this paper, we have studied the back-reaction of long wavelength linear fluctuations on the power spectrum, 
produced by the mode mixing which occurs as a consequence of the non-linearity of the Einstein equations. We 
find that, assuming that the phases of the linear fluctuations are random, there is no phase space enhancement of 
the back-reaction effect. The leading infrared back-reaction contributions are suppressed by neipi compared to the 
contribution of the linear fluctuations, where Ktpi is a measure of the amplitude of the linearized metric fluctuations, 
and e is an inflationary slow-roll parameter. These leading back-reaction terms do not modify the tilt of the power 
spectrum on scales substantially smaller than the Hubble radius. Note that in the case of correlated phases of the 
linear fluctuations, a much larger back-reaction effect is possible. 

We have also seen that the modifications to the Gaussianity of the leading order perturbations are negligible. 
The small size of these modifications goes a long way towards justifying the linear approximation to cosmological 
perturbation theory. However, our results do not exclude the possibility that large amplitude local fluctuations can 
effect the measured fluctuations, as very recently suggested in [3(J (based on the second order formalism developed 
in IH) 0. 

Our work differs from previous work on second order fluctuations in that it emphasizes the fact that, in an acceler- 
ating universe, the phase space of super-Hubble modes is increasing in time. In contrast to what occurs in the case of 
the back-reaction on the homogeneous mode, the back-reaction on the fluctuating modes themselves does not increase 
without limits as a function of time. Compared to previous analyses, our work also gives an easier way to derive the 
leading-order effects of long-wavelength cosmological fluctuations. Our results have been derived in the context of an 
arbitrary slow-roll inflationary model and are thus valid for a wide range of cosmological scenarios. 
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There has been a lot of recent interest in the possibility that the leading gradient terms of long-wavelength modes might 
back- react on local observables in a way that mimics dark energy pfil ] (see also pjfT]'). However, objections to this possibility 
have been raised [3lL 1321 |33J, 134 . l35|| . 
[41] It is also in this approximation that we can justify writing the perturbed metric in the form Q17|l - see e.g. 18] for a 
discussion of this point. 

[42] The possibility that local fluctuations can have a measurable effect on background quantities such as the deceleration 
parameter has recently been suggested in |38l l39|l . 



